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ABSTRACT 

We present new Hubble Space Telescope images of high-velocity H-o and Lyman-o emission in the 
outer debris of SN 1987A. The H-o images are dominated by emission from hydrogen atoms crossing 
the reverse shock. Eor the first time we observe emission from the reverse shock surface well above 
and below the equatorial ring, suggesting a bipolar or conical structure perpendicular to the ring 
plane. Using the Ho imaging, we measure the mass flux of hydrogen atoms crossing the reverse shock 
front, in the velocity intervals (—7,500 < Vohs < —2,800 km s“^) and (1,000 < Vohs < 7,500 km s“^), 

Mh = 1-2 X 10“^ Mq yr“^. We also present the first Lyman-o imaging of the whole remnant and 
new Chandra X-ray observations. Comparing the spatial distribution of the Lyman-o and X-ray 
emission, we observe that the majority of the high-velocity Lyman-a emission originates interior to 
the equatorial ring. The observed Lyman-a/H-a photon ratio, {R{La/Ha)) ~ 17, is significantly 
higher than the theoretically predicted ratio of ~ 5 for neutral atoms crossing the reverse shock front. 

We attribute this excess to Lyman-a emission produced by X-ray heating of the outer debris. The 
spatial orientation of the Lyman-a and X-ray emission suggests that X-ray heating of the outer debris 
is the dominant Lyman-a production mechanism in SN 1987A at this phase in its evolution. 

Subject headings: supernovae: individual (SN 1987A) — shock waves — circumstellar matter 


1. INTRODUCTION 

The reverse shock (RS) in SN 1987A is the sur¬ 
face where the freely expanding debris is suddenly de¬ 
celerated as it encounters circumstellar matter. Hy¬ 
drogen atoms crossing the shock are excited by colli¬ 
sions in the shocked plasma and emit R-a and Lyman- 
a photons having Doppler shifts corresponding to the 
projected velocity of the freely expanding debris im¬ 
mediately i nside the shock . This emission was pre¬ 
dicted by iBorkowski et al.l (Il997f) and subse q uently 
obs erved and interprete d by iSonneborn et al.l (j 19981 ) 
and iMichael et al.l (jl998f ) . 
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The freely expanding supernova debris obeys a “Hub¬ 
ble’s Law”, i.e., the Doppler velocity is given by Vobs = 
zjt , where z is the projected distance along the line 
of sight measured relative to the center of the explosion 
and t is the time since the supernova explosion 27.5 
years in mid-2014). Thus, surfaces of constant Doppler 
shift are planar sections of the debris. The fluxes of H- 
a and Lyman-a photons emitted at the RS are directly 
proportional to the fluence of hydrogen atoms across the 
shock. 

These facts create a unique opportunity to map the 
three-dimensional structure of the RS in SN 1987A and 
provide an empirical framework for studies of the shock 
interaction. This opportunity has been exploited in 
observations with the Space Telescope Imaging Spec¬ 
trograph (STIS). In the STIS observations, the R-a 
emission from the RS was clearly visible near the cir¬ 
cumstellar equatorial ring (ER), where the shock is 
strongest (IMichael et al.l I1998L I2QQ3I: iHeng et ^ I2QQ6I: 
lErance et al.ll2QlQ[) . The ring is inclined ~ 43°, north be¬ 
ing the near side (Panagia et al. 1991; Plait et al. 1995; 
Sugerman et al. 2005). Michael et al. (2003) modeled 
the RS surface, concluding that the observed emission 
was confined to within ^ ± 30° of the ring plane, similar 
to the o pening angle of the ring observed at radio wave¬ 
lengths (jNg et al.ll2013[) . However, the STIS exposures 
were not deep enough to detect R-a emission at high lat¬ 
itudes (i.e., out of the ER plane), where the fluence of 
atoms is much smaller than near the equator. Moreover, 
none of the spectroscopic observations included enough 
STIS slit locations to map the entire shock surface. 

In this Letter, we describe the results of an observing 
campaign to map the high-velocity R-a and Lyman-a 
emission using filters that transmit the blue- and red- 
shifted emission from the debris while suppressing the 
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Fig. 1.— 20 August 2014 STIS G750L spectra with WFC3 fil¬ 
ter curves: dashed cyan lines represent F645N and F665N, while 
the solid red line represents the F656N filter. The F658N filter is 
shifted about 2 nm to the red of F656N and is not shown here. 
Narrow emission lines from the ER are labeled in yellow and the 
broad RS arcs are labeled in orange. The dotted purple line shows 
the position of the equatorial ring. These spectra have an effective 
slit width of 0.5" in the East-West direction of the images shown 
in Figure 2-4. 

bright emission from the ER. In Sections 2 and 3, the 
new HST observations are described; for the first time 
we detect emission from the shock above and below the 
ring plane. In Section 4, we compare H-a, Lyman-a, 
and X-ray observations, arguing that the majority of the 
Lyman-tr flux is driven by X-ray heating of the outer 
ejecta, as proposed in our earlier work (e.g., Larsson et 
al. 2011; France et al. 2011; Fransson et al. 2013). 

2. HST OBSERVATIONS AND IMAGE REDUGTION 

We obtained images of SN 1987A with the WFC3 
camera in four filters to isolate the high-velocity H-a: 
blue high-velocity H-a (“H-a Blue”, F645N), red high- 
velocity H-a (“H-a Red”, F665N), and two narrow-band 
filters centered on shocked emission from the ER and 
hotspots (narrow H-a, F656N and [N II] A 6583 A, 
F658N). The WFC3 images each have exposure times of 
2880 seconds and were acquired on 21 June 2014. Figure 
1 shows the band passes of these filters superposed on 
our 2014 STIS spectrum of SN 1987A (see also Larsson 
et al. - in prep.). 

The “H -a Blue” and “H-a Red” images include a con¬ 
tribution from emission lines originating on the ER (see 
Fransson et al. - in prep, for a discussion of the ER 
emission). The narrow band images are used to subtract 
these contributions to the high velocity material images 
(Figure 2). The F656N image provided the best tem¬ 
plate for the north-dominated high-velocity “H-a Blue” 
image. The F658N filter contains a significant contri¬ 
bution from redshifted H-a on the southern side of the 
remnant, therefore this filter provided the best template 
for the south-dominated high-velocity “H-a Red” image. 
In Figure 2 (right), the ER appears black; however this is 
an artifact of the intentional oversubtraction and stretch 
of the images for display. For the quantitative analysis 
in Section 4, a smaller ER subtraction is applied; the 
emission is nearly continuous across the north and south 
ER region. 

We observed the Lyman-a emission from SN 1987A 


with the ACS/SBC F122M filter for three orbits (9200 
seconds; Figures 3 and 4). The F12 2M filter bandpass 
is shown on an HST-COS spectrum (jFrance et al.ll20lTh 
in Figure 3. In this case, resonant scattering by the neu¬ 
tral ISM of the Milky Way and Large Magellanic Cloud 
acts as a natural blocking filter for low-velocity emission 
from the ER and no additional subtraction is necessary. 
A UV-bright B star (“Star 3”) appears to the south¬ 
east in all imaging bands, providing a common reference 
point. We complement the HST imaging with a recent 
(2014 March 19, day 9885) 0.3 - 8.0 keV image from 
Chandra (Frank et al. - in prep.). 

3. HIGH-VELOGITY H-a IMAGING 

The “H -a Blue” filter samples H-a emission with 
Doppler shifts —7500 < Vobs < — 2800 km s“^, while 
the “H-(a Red” filter samples -1-1000 < Vobs < +7500 km 
s~^ (Figure 2). Regions of interest are labeled A, B, C, 
and D; “-b” and “-r” denote blue and red R-a velocities, 
respectively. Emission from the A-b and A-r regions are 
produced by H atoms passing through the RS near the 
north and south side of the ER, respectively. The bright 
H -a emission features (A-b and A-r) from the RS have 
been partially mapped in 3 dimensions by STIS observa¬ 
tions (e.g., Michael et al. 2003; Heng et al. 2006; France 
et al. 2010). 

In the “H-a Blue” image (Figure 2, top) a bright 
crescent-shaped feature covers the north side of the ER 
(as defined by the [N II] contours), labeled B-b. Given 
the 43° inclination of the ring, the B-b region must orig¬ 
inate from significantly above the ring plane for it to 
be visible in this configuration. This image also shows 
faint streaks that cross the ring extending towards the 
south (C-b). The features labeled A-r, B-r, and C-r are 
analogous to the blue side, but reflected about the ma¬ 
jor axis of the ER. The additional feature labeled D-r is 
probably emission from the internal debris, which is ex- 
pan ding with radial velocities ranging up to ^ 3500 km 
s“^ (jMcCravl 1 19931: [Larsson et al.ll2013[ ). 

Combining the “H-a Blue” and “H-a Red” images, 
we measure the spatial extent of the high-velocity R-a 
regions (positional uncertainties + 0.05"). The high- 
velocity material interior to the ring (A regions) has 
north-south and east-west diameters of 0.94" x 1.31", 
which when including the inclination translates to physi¬ 
cal dimensions of (9.9 x 9.8) x 10^^ cm (assuming d = 50 
kpc). This suggests a symmetric, radial flow with a max¬ 
imum radial velocity of 5.7 x 10^ km s“^ (at 27.5 years 
since the supernova explosion). The B regions on the 
blue and red side have maximum angular distances from 
the center of the remnant of 0.82" and 0.98", respectively. 
These correspond to maximum de-projected distances 
and radial velocities of 8.6 x 10^^ cm and 10.0 x 10^ km 
s“^ for the north/blue and 10.4 x 10^^ cm and 11.9 x 10^ 
km s“^ for the south/red. The C region has a maximum 
observed angular extent of 1.37" from the center of the 
remnant, corresponding to 14.5 x 10^^ cm and 16.7 x 10^ 
km s“^ assuming free expansion. 

3.1. Morphology of the Reverse Shoek Surfaee 

[Fransson et al ] (IMl noted that the large observed 
velocity of the R-a at late times (Vobs ^ H x 10^ km 
s“^ at t > 20 years) suggests that the reverse shock front 
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Fig. 2.— The left column shows the original F645N (“H-a Blue”) and F665N (“H-a Red”) images (21 June 2014), with the blow up 
region outlined in gray. The right column shows an enlarged version after subtracting emission from the equatorial ring. The [Nil] contours 
are shown in magenta to identify the location of the ER. The regions of interest are approximately outlined in dashed white lines. 


has expanded beyond the ring plane. They conclude that 
the RS extent in the polar direction is ^ 15% larger than 
the ER radius and ^ 40% larger than the RS radius in 
the ring plane. We find de-projected velocities of 10 - 
12 X 10^ km s ^ in the B-b and B-r regions, correspond¬ 
ing to RS radii 40 - 70% larger than the 6.1 x 10^^ cm 
ring radius. The C regions extend to 2 - 2.5 times the 
equatorial ring radius. 

To interpret the high-velocity H-a images, we recall 
that surfaces of constant Doppler shift are planar sec¬ 
tions of the freely expanding supernova debris. There¬ 
fore, the pass bands of the “H-a Blue” and “H-a Red” 
filters map to slabs in physical space, as illustrated in 
Figure 5. Portions of the RS surface visible through the 
“H -a Blue” filter (colored blue) reside within the slab 
delineated by the blue vertical lines, while those visible 
through the “H-a Red” filter (colored red) reside within 
the slab delineated by the red vertical lines. The density 
of the outer debris crossing the RS fr ont, and therefor e 
the mass flux, is psn ^ po V' ^ (jLuo et al.l Il994l ). 
The much higher density of the debris crossing the RS 
in the equatorial direction relative to the polar direc¬ 
tion {per/P poi ^ 10^) makes the A region much brighter 
(detectable spectroscopically within ^10 years of the ex¬ 
plosion). The lower density material flowing out of the 
ring plane (B and C regions) is fainter, only now visible 
with the dedicated high-velocity H-a imaging and suffi¬ 
cient time for the polar shock structure to be angularly 


separated from the ring plane. 

While much of the RS surface is not imaged because 
of the discrete velocity surfaces transmitted through the 
filters (shaded gray in Figure 5), the emerging picture 
for the extended H-a emission in the off-plane regions 
suggests a bipolar morphology. This is expected due to 
the high pressure of the equatorial plan e, owing to the 
presence of the ring (jBlondin et al.lll99^ . Assuming the 
outer ejecta pressure is proportional to psn the C po¬ 
lar region pressure is ^ 1/1800 of the equatorial region 
pressure. This general picture resembles the bipolar sur¬ 
faces sugg ested for the pre-supernova interacting wind 
picture bv iBlondin fc LundqvistI (jl993f ). however, recent 
studies of a ‘pre-explosion twin’ to the SN 1987A sys¬ 
tem (SBWl; Smith et al. 2013) suggest that interacting 
winds may not be required to explain the nebular mor¬ 
phology. 

A second possible explanation for the C region emis¬ 
sion is a conical or helical structure (e.g. Smith et 
al. 2007, 2013) connecting the ER with the extended 
outer rings (ORs). In this scenario, 1) high-velocity 
ejecta crossing this structure would be limb-brightened, 
explaining the bright edges in the C-b and C-r emis¬ 
sion, and 2) the southern C-b streaks would be offset 
to the east compared to the northern C-r streaks by 
an amount comparable to the E-W offset of the ORs 
(e.g, Sugerman et al. 2005), which is indeed what is 
observed in Figure 2. In this picture, the C emission 
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Fig. 3. — (top left) 2011 COS spectrum, with F122M filter curve overplotted (dashed red line). At upper right, R(La/Ha) (photon 
flux ratio) in a 0.225^' E-W spatial extraction region (orientation indicated on image below) with normalized X-ray flux in blue circles. 
Emission from the ER prevents a measurement of R{La/Ha) at angles > ± 0.6". R{La/Ha) is elevated towards the regions of brightest 
X-ray emission. The ACS/SBC F122M Lyman-o; image is shown at bottom (June 21 2014). 


would essentially be a new SN 1987A intermediate ring 
system. Deep R-a imaging spectroscopy to map the re¬ 
gion could distinguish between these scenarios and may 
shed light on the progenitor structure, including the pos- 
sibili ty of a binary merger prior to the supernova explo¬ 
sion ([Morris fc PQdsiadlowskill2QQ7D . 

4. HIGH-VELOCITY LYMAN-a IMAGINC 

Figure 3 illustrates the contribution of different far- 
UV emission components to the F122M image. Inter¬ 
stellar H I blocks the center of the Lyman-a line profile 
over the velocity range —15 00 km s“^ < Vobs < +1500 
km s“^ ([France et al.l [2QT11) . and so excludes emission 
from the shocked ring and hotspots, which have Doppler 
shifts A V < 300 km s“^ (Pun et al 2002). Based on 
the HST-COS spectrum, we estimate that 90 - 95% of 
the total emission in the F122M image is contributed by 
high-velocity Lyman-a. The remaining 5 - 10% is mostly 
N V A 1240 A emission from the ER and high-velocity 
ions crossing the RS front. 


Separation of the peak Lyman-a and ER emission is 
complicated on the bright northern side by the 43° pro¬ 
jection on the plane of the sky, however, there is separa¬ 
tion between the Lyman-a and the R-a ring profiles in 
the East-West direction. The East-West ring diameters 
are 1.44" (± 0.05"), 1.62" (± 0.05"), and 1.63" (± 0.10") 
for Lyman-a, R-a ER (E656N), and X-ray images, re¬ 
spectively. The Lyman-Of image is slightly offset with re¬ 
spect to the ER, with the largest separation (1.35 x 10^^ 
cm) at the western edge (Eigure 4, bottom). While the 
X-ray resolution is not as high as the HST imaging of 
the ER, there is a clear offset between the X-ray and 
Lyman-Of profile peaks. Unlike the X-ray and Lyman-a 
images, wh ich are brightest in t he NW quadrant, the ra¬ 
dio image ([Zanardo et al.l[2014f ) is brightest in the East. 

4.1. The Lyman-a/H-a Ratio: X-ray Heating 

Neutral Hydrogen Mass Flux and the Lyman-a/H-a 
Ratio - If all of the high-velocity hydrogen emission 
originates from collisional excitation as neutral hydro- 
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Fig. 4.— ACS/SBC F122M Lyman-o: image with 19 March 2014 
Chandra (0.3 - 8.0 keV) contours overlaid in blue (top). Below, 
spatial profiles of the western edge of the ER region show that the 
Lyman-o; emission arises interior to the equatorial ring as traced 
by X-ray (blue), [O III] (green), and low-velocity H-o emission 
(black). 


gen atoms cross the RS front, the total flux is propor¬ 
tional to the fluence of hydrogen atoms across the shock 
boundary. Integrating the individual band images, we 
calculate the total number of atoms in a given velocity 
range. We take the total (Galactic + LMC) reddening 
towards SN 1987A as a standard ISM curve with Ry = 

3.1 and E{B — V) = 0.19 (see France et al. 2011 for 
a discussion). The total rate of hydrogen atoms crossing 
the RS front is approximately: 

Nh = X F(AA)-^AA x P{H)-^ x Cjsm (1) 
he 

where d ^ 50 kpc (^ 10% uncertainty; Panagia et al. 
1991), F{AX) is the image-integrated flux density in a 
given band (^ 5% uncertainty), AA is the effective band¬ 
pass of the filter, P{H) is the number of line photons 
emitted per atom, and Cjsm is the interstellar attenu¬ 
ation corr ection factor. P(H) ^ 0.2 for H-a and 1 for 
Lyman-(a (|Michael et al.|[2QQ^ iHeng fc McCravll2007D . 

For the “H-a Blue” image, AA 100 A and Cjsm 
= 1.57 ± 0.08, giving a total hydrogen fluence of 

4.1 (± 0.7) X 10^^ atoms s“^. For the “H-a Red” 



Fig. 5. — Schematic representation of the underlying structure 
in the H-o; imaging. The diagonal lines connecting the filled circles 
represent the circumstellar rings. The red (blue) vertical lines de¬ 
lineate the slabs in which H-a emission from the reverse shock is 
visible through the “H-a Red” ( “H-a Blue”) filter. The A, B, C, D 
labels refer to the image locations in Figure 2; the “b” and “r” refer 
to the “H-a Blue” and “H-a Red” images, respectively. Regions 
colored blue and red appear in Figure 2, while regions colored gray 
are invisible to the H-a imaging. 


image, AA « 120 A and Cjsm = 1.53 ± 0.08, giving 
a total hydrogen fluence of 3.3 (± 0.5) x 10^^ atoms 
s“^. Taken together, we estimate the total mass flux of 
high-velocity hydrogen atoms crossing the RS is Mh = 
7.4 (± 1.2) X 10^2 g s-i (1.2 X 10-^ Mo yr-^). 

The analogous calculation using the F122M Lyman- 
a image has AA « 100 A and Cjsm = 6-7 x 1.63, 
where 6.7 (± 1.7) is the correction for dust attenuation 
and 1.63 (± 0.08) is the empirically determined correc¬ 
tion for the interstellar Lyman-a line core attenuation. 
This calculation shows a Lyman-a-to-H-a photon ratio, 
R{La/Ha), of ^ 17 ± 6, a factor of 3.5 larger than 
attributable to excitation in the RS alone. Therefore, 
an additional Lyman-a emission mechanism is present. 
Note that the reddening correction does not account for 
differential extinction within the supernova debris, so the 
intrinsic R{La/Ha) may be greater. Lyman-a excess 
was observed by previous authors using spectroscopic ob¬ 
servations o f small spatial regions of the SN 1987 A inner 
ring region (|Heng et al.ll2QQ6HFrance et al.llMTol) . 

Combining the “H-a Blue” and “H-a Red” images with 
the Lyman-(a map, we And the baseline R{La/Ha) in¬ 
terior to the ER is 8 - 10, at the center and southeast¬ 
ern sides of the inner remnant. R{La/Ha) reaches a 
maximum on the western and northwestern side of the 
inner remnant where the Chandra emission is bright¬ 
est (Figure 3). The maximum R{La/Ha) is ~ 35 near 
PA = 270°, just interior to the brightest X-ray emission 
from the ER. R{La/Ha) values ^ 20 are found towards 
PA ^ 60°, the location of the northeast X-ray maximum. 

The Role of X-rays - The observed offset between 
the X-ray and Lyman-a emitting regions is consistent 
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with the expected separation between the ER and the 
outer extent of the RS. Following the observation that 
the dominant energy input in the unshocked supernova 
debris has evolved from internal heating by radioac¬ 
tive dec ay to external heati n g from X-rays produc ed at 
the ER (jLarsson et al.ll2Ql'Tl) . iFransson et alJ (j2Q13[ ) pre¬ 
sented calculations of the X-ray energy deposition into 
the ejecta. The heating of the outer debris is dominated 
by X-rays with E < 0.5 keV while higher energy X-rays 
penetrate deeper to heat the interior. X-rays propagat¬ 
ing into the outer debris produce fast photoelectrons, 
which deposit most of their kinetic energy as heat, pro¬ 
vided thn^_thn_d£teisji^ fractional ionization n^juH > 
0.03 (|Xn fc McCravIfToMI ). Radiative cooling by thermal 
excitation of Lyman-a limits the debris temperature to 
T < 10^ K, at which temperature thermal excitation 

of H-Of is negligible (this mechanism, and not recombi¬ 
nation, also accounts for the H I two-photon spectrum 
observed in the UV, France et al. 2011). Therefore, 
we expect the Lyman-a/H-a ratio to be at a maximum 
nearest the source of the X-ray heating. The spatial 
stratification of the Lyman-a and X-ray emission and 
the enhanced Lyman-a/H-a ratio support the scenario 
where X-ray heating of the outer debris is responsible 
for the majority of the observed Lyman-a emission from 
SN 1987A at present. 

The total Lyman-a luminosity in the F122M band is 
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^hyoL = 2.5 (± 0.7) X 10^^ erg ^ 70% of which is in 
excess of the expected RS emission. The total flux of the 
X-ray ring is 1.04 x 10“^^ erg cm“^ s“^, or, Lx (0.3 - 8.0 
keV) = 3.1 (± 0.3) X 10^^ erg s“^. However, most of the 
heating radiation is in the soft X-ray/EUV band (0.01 - 
0.5 keV) which is not observed at Earth because of in¬ 
terstellar neutral hydrogen and dust attenuation towards 
SN 1987A. Scaling the Zhekov et al. (2006) model X-ray 
spectrum to the observed day 9885 0.3 - 8.0 keV flux, 
the total 0.01 - 8.0 keV luminosity is Lx (0.01 - 8.0 keV) 
= 5.6 X 10^^ erg s“^. Assuming that the only source of 
EUV/soft X-rays is that observed by Chandra and that 
the neutral hydrogen in the outer ejecta intercepts half of 
the soft X-rays emitted by the ring, the Lyman-a heat¬ 
ing efficiency must be ^ 60% to reproduce the F122M 
image of SN 1987A. However, slower shocks that con¬ 
tribute more to the optical/UV emission of the ER than 
to the X-rays wi ll also be a source of EUV irradiance of 
the outer debris (jEransson et al.l[2013l ) , meaning that the 
actual soft X-ray heating efficiency is likely substantially 
lower. 

The HST observations presented here were acquired as 
part of the Cycle SAINTS program (13401 and 13405). 


The HST observations presented here were acquired as 
part of the Cycle SAINTS program (13401 and 13405). 
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